对流扩散方程的QC3无网格法

对流扩散方程作为偏微分运动方程的分支,在流体力学、气体动力学等领域有着重要应用.为解决对流扩散方程难以通过解析法得到解析解的难题,采用二阶一致3点积分(Quadratically Consistent 3-Point Integration,简称QC3)提高无网格法的计算效率,通过对积分点上形函数导数的修正,改善无网格法的精度和收敛性.本文将QC3无网格法拓展到对流扩散方程问题中,时域离散采用广义特征线Galerkin法,空间离散采用QC3法.数值结果表明,应用QC3无网格法得到的对流扩散问题数值解与解析解十分接近,验证了QC3无网格法解决对流扩散问题的可行性.     

移动加热器法生长碲锌镉晶体的组分输运与界面形貌研究

移动加热器法(THM)生长碲锌镉晶体时,界面稳定性对晶体生长的质量有很大影响。本文基于多物理场有限元仿真软件Comsol建立了THM生长碲锌镉晶体的数值模拟模型,讨论了Te边界层与组分过冷区之间的关系,对不同生长阶段的物理场、Te边界层与组分过冷区进行仿真研究,最后讨论了微重力对物理场分布的影响,并对比了微重力与正常重力下的生长界面形貌。模拟结果表明,Te边界层与组分过冷区的分布趋势是一致的,在不同生长阶段,流场中次生涡旋的位置会发生移动,从而导致生长界面的形貌随着生长的进行发生变化,同时微重力条件下形成的生长界面形貌最有利于单晶生长。因此,在晶体生长的中前期,对次生涡旋位置的控制和对组分过冷的削弱,是THM生长高质量晶体的有效方案。     

从化学视角剖析烫发、染发原理及顺序*

从头发的结构和组成出发分析头发形状和颜色改变的可能性,剖析烫发、染发的化学原理,阐释先烫发后染发的本质原因,从化学视角辨证地看待烫发、染发的利弊。     

Chemical bonding in representative astrophysically relevant neutral,cation, and anion HCnH chains

Most existing studies assign a polyynic and cumulenic character of chemical bonding in carbon-based chains relying on values of the bond lengths. Building on our recent work, in this paper we add further evidence on the limitations of such an analysis and demonstrate the significant insight gained via natural bond analysis. Presently reported results include atomic charges, natural bond order and valence indices obtained from ab initio computations for representative members of the astrophysically relevant neutral and charged HC_2k/2k+1H chain family. They unravel a series of counter-intuitive aspects and/or help naive intuition in properly understanding microscopic processes, e.g., electron removal from or electron attachment to a neutral chain. Demonstrating that the Wiberg indices adequately quantify the chemical bonding structure of the HC_2k/2k+1H chains—while the often heavily advertised Mayer indices do not—represents an important message conveyed by the present study.     

Effect of wall slip on the viscoelastic particle ordering in a microfluidic channel

The formation of a line of equally spaced particles at the centerline of a microchannel, referred as “particle ordering,” is desired in several microfluidic applications. Recent experiments and simulations highlighted the capability of viscoelastic fluids to form a row of particles characterized by a preferential spacing. When dealing with non-Newtonian fluids in microfluidics, the adherence condition of the liquid at the channel wall may be violated and the liquid can slip over the surface, possibly affecting the ordering efficiency. In this work, we investigate the effect of wall slip on the ordering of particles suspended in a viscoelastic liquid by numerical simulations. The dynamics of a triplet of particles in an infinite cylindrical channel is first addressed by solving the fluid and particle governing equations. The relative velocities computed for the three-particle system are used to predict the dynamics of a train of particles flowing in a long microchannel. The distributions of the interparticle spacing evaluated at different slip coefficients, linear particle concentrations, and distances from the channel inlet show that wall slip slows down the self-assembly mechanism. For strong slipping surfaces, no significant change of the initial microstructure is observed at low particle concentrations, whereas strings of particles in contact form at higher concentrations. The detrimental effect of wall slip on viscoelastic ordering suggests care when designing microdevices, especially in case of hydrophobic surfaces that may enhance the slipping phenomenon.     

Subset representations and eigenvalues of the universal intertwining matrix

We solve a combinatorial question concerning eigenvalues of the universal intertwining endomorphism of a subset representation.     

Shedding vortex simulation method based on viscous compensation technology research

Owing to the influence of the viscosity of the flow field, the strength of the shedding vortex decreases gradually in the process of backward propagation. Large-scale vortexes constantly break up, forming smaller vortexes. In engineering, when numerical simulation of vortex evolution process is carried out, a large grid is needed to be arranged in the area of outflow field far from the boundary layer in order to ensure the calculation efficiency. As a result, small scale vortexes at the far end of the flow field cannot be captured by the sparse grid in this region, resulting in the dissipation or even disappearance of vortexes. In this paper, the effect of grid scale is quantified and compared with the viscous effect through theoretical derivation. The theoretical relationship between the mesh viscosity and the original viscosity of the flow field is established, and the viscosity term in the turbulence model is modified. This method proves to be able to effectively improve the intensity of small-scale shedding vortexes at the far end of the flow field under the condition of sparse grid. The error between the simulation results and the results obtained by using fine mesh is greatly reduced, the calculation time is shortened, and the high-precision and efficient simulation of the flow field is realized.